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ABSTRACT 

Doubly-charged Higgs bosons (A^^/A++) appear in several 
extensions to the Standard Model and can be relatively light. We 
review the theoretical motivation for these states and present a 
study of the discovery reach in future runs of the Fermilab Teva- 
tron for pair-produced doubly-charged Higgs bosons decaying 
to like-sign lepton pairs. We also comment on the discovery 
potential at other future colliders. 

I. Introduction 

Doubly-charged Higgs bosons (A^^) appear in exotic Higgs 
representations such as found in left-right symmetric models. 
The current experimental bound is m^-- > 45 GeV from 
a search for Z" A — A++ at LEP. 

At the Tevatron, the two production mechanisms with poten- 
tially large cross section are pair production, pp — > -f/Z'^X 
A A++X or single production via WW fusion, pp 
W~W~X A~~X. However, existing phenomenologi- 
cal and theoretical constraints are only easily satisfied if the 
W~W^ A^^ coupling is vanishing (or very small). There- 
fore, in this analysis we will consider the discovery reach for 
detecting A^^ A++ pair production at the Tevatron. 

In many models, it is possible for the A to couple to like- 
sign lepton pairs, £^£^ . If the W^W^ — > A^^ coupling is 
vanishing, it is then very likely that the doubly-charged Higgs 
will decay to £^£^ via the lepton-number-violating coupling. 
We will therefore concentrate upon A"~ — > e~e~, A^^ 
and A t^t^ . 

Alternatively, if the A — -> £-£- and A — W^W~ 
couplings are both vanishing or very small, then the A can 
have a sufficiently long lifetime that it will decay outside the 
detector. Identification of the A A++ pair via the associated 
dE/dx distributions in the tracking chamber would then be pos- 
sible. 
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11. Theoretical Motivation 

Doubly-charged scalar particles abound in exotic Higgs rep- 
resentations and appear in many models ||, ^. For example, 
a Higgs doublet representation with Y = —3 contains a doubly- 
charged A and a singly-charged A^. If part of a multiplet 
with a neutral member, a A^^ would immediately signal the 
presence of a Higgs representation with total isospin T = 1 or 
higher Most popular are the complex Y = —2 triplet Higgs 
representations, such as those required in left-right symmetric 
models, that contain a A , a A^ and a A°. 

In assessing the attractiveness of a Higgs sector model con- 
taining a A^^ many constraints need to be considered. For 
triplet and higher representations containing a neutral member, 
limits on the latter's vacuum expectation value (vev) required 
for p = rn^ /[cos^ 6 wrn^z] = 1 at tree-level are generally 
severe. (The first single representation beyond T = 1/2 for 
which p = 1 regardless of the vev is T = 3, F = —4, whose 
T3 = member is doubly-charged.) Models with T = 1 and 
T = 2 can have p = 1 at tree-level by combining represen- 
tations. However, such models generally require fine-tuning in 
order to preserve p = 1 at one-loop. The simplest way to avoid 
all p problems is to either consider representations that simply 
do not have a neutral member (for example, a F = — 3 doublet 
or a y = —4 triplet representation), or else models in which 
the vev of the neutral member is precisely zero. We will only 
consider models of this type in what follows. 

Further constraints on Higgs representations arise if we re- 
quire unification of the coupling constants without intermedi- 
ate scale physics. In the Standard Model, unification is possi- 
ble for a relatively simple Higgs sector that includes a single 
|y| = 2 triplet in combination with either one or two |y| = 1 
doublets (the preferred number of doublets depends upon the 
precise value of as{mz))- In the case of the minimal super- 
symmetric extension of the Standard Model, precise unification 
requires exactly two doublet Higgs representations (plus pos- 
sible singlet representations); any extra doublet representations 
(including ones with a doubly-charged boson) or any number 
of triplet or higher representations would destroy unification. 
However, by going beyond the minimal model and including 
appropriate intermediate-scale physics, supersymmetric models 



(in particular, supersymmetric left-right symmetric models [Q]) 
with triplet and higher representations can be made consistent 
with unification. 

In short, the popular two-doublet MS SM need not be nature's 
choice. We should be on the look-out for signatures of exotic 
Higgs representations, the clearest of which would be the exis- 
tence of a doubly-charged Higgs boson. Thus, it is important to 
consider how to search for and study such a particle. 

The phenomenology of the A^^ derives from its couplings. 
Tri-Unear couplings of the type W^W^ A^^ are not 
present in the absence of an enabling non-zero vev for the neu- 
tral member (if present) of the representation, and q'qA cou- 
pUngs are obviously absent. There are always couplings of the 
form 7 ^ A A++. In addition, and of particular inter- 
est, there is the possibility of lepton-number-violating £^£^ 
A couplings in some models. For (5 = 73 + ^ = ^2 the 
allowed cases are: 



^L^R 



A- 

A- 
A- 



■(T = o,r3 = o,y = -4), 

-(T=i,T3 = -i,r=-3). 

-(r = i,r3 = -i,r = -2). 



(1) 



Note that the above cases do not include the T — 2>,Y ^ —4 
representation that yields p — \, nor the T — 1,Y = —A triplet 
with no neutral member, but do include the T — 1/2, F = —3 
doublet representation with no neutral member, and the popular 
T ~ \,Y — —2 triplet representation. In left-right symmet- 
ric models there is a 'right-handed' and a ieft-handed' Higgs 
triplet, both with |y| — 2. Our analysis applies to the left- 
handed triplet (whose neutral member must have a very small 
vev to preserve p = 1); the phenomenology of the right-handed 
triplet is completely different. 

In the case of a |y| = 2 triplet representation (to which we 
now specialize) the lepton-number-violating coupling to (left- 
handed) leptons is specified by the Lagrangian form: 



Cy = ihijiiJj^CT2AipjL + h.c. 



(2) 



where i, j — e, p,T are generation indices, the ^/I's are the two- 
component left-handed lepton fields {iI^il = {i^i,^^ and A 
is the 2x2 matrix of Higgs fields: 



A = 



A-/V2 
AO 



A— 
-A-/V2 



(3) 



Limits on the hij coupling strengths come from many 
sources. Experiments that place limits on the hij by virtue of the 
A £~£^ couplings include Bhabha scattering, {g — 2)^, 

muonium-antimuonium conversion, and /i^ e^e^e+. These 
limits [|], ^ suggest small off-diagonal couplings (as assumed 
in our analysis). Writing 



naturally). Single production of A , A++ (A ) is possible 
in e+e" (ep) collisions at LEP2 (HERA) via diagrams involv- 
ing the A^^ er or A++ e+e+ couplings. If 
saturates its upper limit, then LEP2 and HERA will probe up to 
TO^-- ^ 150 GcV [H^. However, it is likely that Cgg is much 
smaller than its current bound and that these sources of single 
production will be negligible. 

Thus, we focus on 7*, — > A^^ A++ pair production, the 
cross section for which is determined entirely by the quantum 
numbers of the A^^. For a general spin-0 boson B, with weak 
isospin T3 and charge Q, and a fermion /, with ^3 and q, the 
// BB pak-production cross section is: 



xwVw 



-P: 



zz- 



(5) 



where s is the // center of mass energy squared, /3 = 

^ xw,A = T:i-xwQ, 



\/\~ 4m'g/s, x\Y — sin^ 9w, Vw = 1^ 
flL = <3 - xwq, an = -xwq, -P77 



t2 \2 



^'z^'z\ 



', an = -Xwq, P-jf = s , Pzz = [{s - 
and Pjz = {s — Tn?z)Pzz / s. We will con- 
sider a A with Ts = -1, Q = -2. An extra factor of 1/3 
is required for color averaging in qq annihilation in pp or pp 
collisions. 




sqrt(s) = 1 4 TeV 



_Tevatron 
sqrt(s) = 2 TeV 



A mass (C^eV) 



Figure 1 : A++ / A pair production cross section as a function 
of A^^ mass for both the Tevatron and the LHC. 
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the limits imply c^e ^ ^.^^ 

Regarding production mechanisms, the fusion process [|], ^ 
0], W-iy- ^ A- , is absent since the required tri-linear cou- 
pling is zero if the vev of the neutral member (if there is one) 
of the Higgs representation is zero (as we assume so that p = 1 



In e+e^ A^^A++, kinematic reach is limited to 

rri/x-- < y/s/2,i.e. no more than about 230 — 240 GcV at a fu- 
ture = 500 GeV NLC. We will find that the discovery reach 
at the Tevatron can cover much, if not all, of this range, depend- 
ing upon the dominant A~~ decay mode. The mass reach for 
pair production at a pp collider increases rapidly with machine 
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energy. Figure [l| shows the A A++ pair production cross 
section for both the Tevatron (at ^ = 2 TeV) and the LHC. At 
the Tevatron, a^''" - 0.9(0.24) fb at m^— = 250(300) GeV. 
With total accumulated luminosity of 30 fb^^ (as planned for 
the TeV33 upgrade) there would be about 27(7) A A"''+ 
events. The marginality of the latter number makes it clear that 
m^-- < 300 GeV will be the ultimate mass reach possible at 
the Tevatron. 

Decays of a A are generally quite exotic ^. For ^ 
A — > W~W~ coupling, the only two-body decays that 
might be important are A — > A^M^^, A^^ A^A~ 
and, if the lepton coupling is present, A . Typically, 

the A^^ and A^ have similar masses, in which case A 
A~A~ is likely to be disallowed. Thus, we will focus on the 
A~M^" and i~ final states. For a T = 1, y = -2 triplet we 
find r " 



^ A-- 
^ A-- 



167r 
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(1.3 GeV) 



100 GoV. 



ImA- ~(0.4GeV)(^)( 



, 100 GcV; 



(6) 

where (3 is the usual phase space suppression factor, and 
we used Eq. (Q). For example [^, if toa-- = 360 GeV, 
mA- = 250 GeV we find r(A— ^ A^V/-) - 2 GeV and 
r(A— ^ £-£-) = 19 GeV (^). If any c« is near IQ-^ 
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limits on Ct-t, the r r 



is likely. Since there are currently no 
channel could easily have the largest 



partial width and be the dominant decay of the A^^. On the 
other hand, if all the cu are very small then the A^W^ mode 
is quite likely to be dominant if it is kinematically allowed. The 
implications for detection of A^^A++ pairs will now be dis- 
cussed. 

in. Simulation and Reconstruction 

The signal and backgrounds are simulated with the PYTHIA 
Monte Carlo, which has been modified to allow the process: 



PP 



Z°/"fX ^ A—A++X, 



(7) 



with the A^^ then forced to decay to like-sign lepton pairs. 

The events are then fed to a CDF detector simulation which 
includes the geometry of the Run I CDF detector. For the Main 
Injector runs of the Tevatron, the CDF and DO detectors will 
both be upgraded to handle higher instantaneous luminosity. 
In addition, the acceptances of the upgraded detectors will im- 
prove. This simulation includes muon coverage for |?7| < 1, 
which will be improved to \rj\ < 1.5 for Run II. This results in 
approximately a 20% improvement in acceptance for this pro- 
cess. 

Events are passed through the normal CDF event reconstruc- 
tion package. Muon candidates must have tracks in both the 
central tracking and the muon chambers, electron candidates 
must have a track and an isolated electromagnetic calorimeter 
cluster. The lepton momenta are determined from the central 
tracking chamber and, if fiducial, the silicon microvertex detec- 
tor For tracks which do not pass through the microvertex detec- 
tor, the fit is performed assuming that the track originated from 



the interaction point. This so-called "beam-constraint" signifi- 
cantly improves momentum and, hence, mass resolution. 



IV. A 



e e 



For the case where the A decays to like-sign leptons (exclud- 
ing taus), the signature is a spectacular 4e or 4/i final state. Here 
we will focus upon the 4/i final state. Backgrounds are very sim- 
ilar for the two channels, although the discovery reach will be 
slightly higher in the electron channel due to better mass reso- 
lution and larger electromagnetic calorimeter coverage. 

The dominant backgrounds in the 4/i mode (accepting at 
least 2 same-sign /i's as described below) arise from elec- 
troweak processes where real high-pT muons are created from 
W or Z decays along with either fake muons or muons from 
heavy flavor decay. The backgrounds are diboson production 
(ZZ 4/i, WZ^ 3/i + I/, WW ^ 2/i + 2iyy, tt production 
(ti — *■ n'^iyb ii~Vhy, and boson plus jets {W + jets, Z + jets), 
where W ^iv, Z /^^/i^ and the jets produce real or fake 
muons. We use the measured cross sections for tt, W + jets 
and Z + jets [0 [u], [l|] and the calculated cross sections for 
W Z and ZZ production [Q. The PDF world average branch- 
ing ratios are used for Z /i^A*~ (0.03367) and W ii^Vi,. 
(0.104) 
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Figure 2: Background contributions to the same sign mass 
plot after all cuts. As can be seen in the inset, the dominant 
background above 100 GeV is from diboson production. The 
Njct < 1 jet requirement removes most of the background. 



With a small event rate and small backgrounds, it is desirable 
to keep the selection criteria as loose as possible. To simulate 
the trigger, we require one central muon with > 12 GeV/c. 
Next, at least two same-sign muons are required. This is suf- 
ficient to easily isolate a signal in the low mass region m^^ < 
100 GcV, but the background from dimuon top decay becomes 
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Figure 3: The effective cross section (true cross section x ef- 
ficiency X acceptance) for signal (unshaded) and background 
(shaded) for A^^ masses of 100, 200 and 300 GeV respec- 
tively, after all cuts have been applied. 



an issue at higher masses. We expect very little additional activ- 
ity in A pair production other than the energy recoiling against 
the virtual Z/j. For dimuon top decay backgrounds, there is 
additional jet activity from the two b decays. The third muon 
is supplied by one of these b decays. The background in the 
high mass region from top decays can be greatly suppressed by 
requiring that the event has no more than one jet seen in the 
calorimeter (\ri\ < 2.5) with more than 7.5 GeV of transverse 
energy (Et). 

Figure || shows the same-sign muon mass contributions from 
each of the backgrounds listed above after all cuts. Above 
100 GeV in same-sign mass, the dominant background is di- 
boson production. 

Figure ^ shows the same-sign dimuon invariant mass distri- 
butions for both signal and background for three different A 
masses: 100, 200 and 300 GeV after all cuts have been ap- 
plied. The signal-to-background ratio remains high at m^-- = 
300 GeV, although the dimuon mass resolution is worsening. 
The dimuon mass resolution is dominated by the pT resolution 
of the detector, which worsens at higher momenta. For the case 
A e^e^ the dielectron mass resolution does not degrade 

as rapidly due to the energy resolution of the electromagnetic 
calorimeter In the case of a high mass search for A~^ decay- 
ing to muon or electron pairs, the technique would be a count- 
ing experiment, looking for an excess of high mass pairs over 
the small background. 

The simulated data shown in Figures ^ and |] represent the 
response of the Run I CDF detector. The product of the effi- 
ciency and acceptance for a signal event to produce at least one 
entry in the same-sign mass plot depends upon rri/^--, but is 



typically 50%. If we assume that the efficiency will be the same 
for the Run II detector, scaling the acceptance to the improved 
muon coverage {\ri\ < 1.5) brings this number up to 60%, an 
improvement factor in the accepted signal of 1.2. The corre- 
sponding scaling of the acceptance for the background is found 
to be ^ 2 over the entire mass region above 50 GcV. 

In addition to the significance of a bump in the same-sign 
mass distribution, there is additional information in the number 
of high-pT muons in the event. With the Run I CDF detector 
used in this simulation, approximately 20% of the signal events 
have four found muons, yielding two entries in the mass distri- 
bution, while ^ 6% of the background events have more than 
three muons such that both same-sign pairs have mass greater 
than 50 GeV. For the Run II detector, these numbers go up to 
approximately 40% for signal and 11% for the background. 
The probability that a background event would have four muons 
and both same-sign combinations near one-another in mass is 
exceedingly small. We therefore conclude that production of 
enough events so that two events are measured to have four 
muons (in addition to the other same-sign dimuon mass entries 
from 2 and 3 muon events) will be more than adequate to es- 
tablish a signal for the A^^. As an example using the numbers 
above: if 10 events are produced, 6 events would produce at 
least one same-sign mass pair Of those 6 events, 2 (from 2.4) 
would have four found muons and 4 (from 3.6) would have 2 
or 3 found muons, yielding 8 entries in the same-sign mass plot 
on a background of approximately one same-sign dimuon mass 
pair and zero four muon events. Tri-muon events offer little ad- 
ditional evidence for A production, since background events 
often have two real, opposite-signed leptons in addition to one 
lepton from either a fake or heavy flavor decay. 

Using the criteria that 10 pair-produced events would lead 
to an unambiguous discovery of the A , we conclude that 
a reach of approximately 200, 250, 300 GeV in the mass of 
the A^^ could be achieved in 2, 10, 30 fb~^ of Tevatron run- 
ning, respectively, for the cases where A e~e^ and 
A fi^ /.i^ . 

V. A— T-T- 

Unlike the electron and muon channels, reconstructing an in- 
variant mass in the tau channel is problematic because of the 
neutrinos involved in their decay. We therefore use a counting 
method to estimate the reach of a doubly charged Higgs search 
for Run II. 

Tau lepton identification is not trivial at a hadron collider. 
Identification efficiencies are much lower than for electrons or 
muons (^ 50%) and fake rates from QCD jets are significant 
(^ 0.5%). Nonetheless, searching the tau lepton channel is 
worthwhile because the doubly charged Higgs may preferen- 
tially couple to the taus, and the tau lepton offers the possibility 
of measuring the spin of its parent. 

Selection of tau lepton candidates which decay into hadrons 
is detailed in [p^. The algorithm begins by looking for jets 
in the calorimeter The tau candidate must have one or three 
charged particles in a 10° cone about the jet axis and no addi- 
tional charged particles in a cone of 30°. In addition, the tau 
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Table I: Expected A tt events passing all cuts in 10 fb 

Ma-- (GeV) events 

50 19. ±4. 

100 8.8 ±0.6 

150 3.07±0.20 

200 0.72±0.11 

250 0.23±0.03 



candidate must have the correct total charge (±1) and have a 
mass consistent with a tau lepton. At a minimum, the cluster 
must have Et > 10 GeV and the largest px of an associated 
charged particle must exceed 10 GeV/c. Finally, the tau candi- 
date cannot be consistent with an electron. 

Currently, a fiducial cut of \ri\ < 1 is required to maintain 
good charged particle tracking efficiency. The tracking cover- 
age for the Run II detector will be significantly larger and tracks 
will be measured with improved resolution. The corresponding 
gains in tau lepton acceptance have not been included in the 
results below. 

Three types of triggers are considered for this search: inclu- 
sive electron and muon triggers for the case where at least one 
of the taus has decayed leptonically, an trigger which relies 
on the neutrinos in the tau decay, and a dedicated tau lepton 
trigger which identifies tau leptons from tracking or calorimetry 
information early in the triggering system. These triggers corre- 
spond to triggers used in Run I and allow the background from 
fake taus to be estimated from data. For the estimates below the 
logical OR of the missing Et and the inclusive lepton triggers 
has been used. The threshold on the trigger is 35 GeV and 
the threshold on the inclusive lepton triggers is Et > 20 GeV. 
The primary tau lepton in an event must have Et > 20 GeV. 

Pair production of doubly charged Higgs produces events 
with four taus in the final state. Approximately 60% of these 
events contain three or more taus which decay into hadrons. 
To reduce the background from fake taus to a reasonable level, 
a selected event must have at least three identified hadronic 
tau lepton candidates and an additional jet, electron, muon, or 
hadronic tau candidate. The expected number of signal events in 
a 10 fb^^ sample for various masses is shown in Table ^ 

The backgrounds are expected to come from two sources — 
processes which produce real taus and processes which produce 
jets which fluctuate to imitate tau leptons. The expected num- 
ber of events containing multiple fake tau leptons was estimated 
from the Run I data sample. No events pass the selection re- 
quirements. 

Top and diboson production are expected to be the largest 
source of background events with real tau leptons. Both of these 
were estimated from Monte Carlo. Again, no events passed the 
strict topology cuts. Given the large data samples expected in 
Run II, both of these backgrounds should be measured rather 
precisely. 

Given that no background events pass the selection require- 
ment, it is difficult to quantitatively define the number of events 
necessary to claim a discovery. Nonetheless, the background is 



likely to be quite small, so a handful of events should be consid- 
ered significant. Arbitrarily taking five events as the standard, a 
search for doubly charged Higgs in the tau channel would have 
a reach of approximately 130 and 180 GeV in samples of 10 
and 30 fb^^, respectively. The cuts, trigger, etc. have not been 
optimized, so these should be considered conservative estimates 
of the Run II reach. 

VI. Conclusions 

Although currently out of favor because of the success of the 
minimal supersymmetric model, there are well-motivated mod- 
els containing triplet and other Higgs representations which in- 
clude a A^^ Higgs boson that has small (most naturally zero) 
W~W~ coupling but possibly non-zero £^£^ coupling. It is 
then very possible that B{A^^ i^i") ^ 1 for £ — e, ^, 
or (most probably?) r. We have demonstrated that detection 
of the A^^ at the Tevatron (operating at y/s = 2TeV with 
i = 30 fb^^) will then be possible for m^-- up to 300 GeV 
for £ = e or fj. and 180 GeV for £ — t. We can estimate from 
Fig. |l]the corresponding limits at the LHC by requiring the same 
raw number of events before cuts and efficiencies as needed at 
the Tevatron — ^10 for £ — e,fi and ^ 300 for £ — t — 
yielding m^-- discovery up to roughly 925 GeV (1.1 TeV) for 
£ ^ e,n and 475 GeV (600 GeV) for £ ^ t, assuming total 
integrated luminosity of L = 100 fb"^ (L = 300 fb"^). For 
£ = e, /i, the reach of the LHC detectors will likely be even 
greater than this, due to the improved lepton acceptance and 
resolution anticipated over the current generation of hadron col- 
hder detectors. For £ = t, this simple extrapolation may not 
account for a different signal-to-background ratio in t selection 
at the LHC. A full study is necessary to evaluate this. 

As detailed in ijl^], if a A is found then e~e" and /i^/i^ 
colliders capable of high luminosity at y/s — m^-- will be- 
come a priority in order to actually determine the cu'&- In- 
deed, observation of A A++ pair production in only a sin- 
gle A £~£~ channel provides no information on c«. (Of 
course, if more than one ££ channel is seen, ratios of the cu's 
could be obtained.) Only if the A — > A'W" decay channel 
[for which the partial width can be computed and compared to 
the £^£~ partial width via Eq. (^] is also seen, can one get an 
estimate of the c« magnitude(s). In contrast, an e~e~ {p,~ 
collider would provide a direct measurement of Cee {c^^i). For a 
more detailed discussion see [p^. This illustrates an important 
complementarity between the NLC and hadron colliders. Dis- 
covery of a A^^ prior to the construction and operation of the 
e+e^, e^e^ collider NLC complex would be very important in 
determining the energy range over which good luminosity and 
good energy resolution for e^e^ collisions should be a priority. 

REFERENCES 

[1] P.D. Acton et al., Phys. Lett. 295B, 347 (1992). 

[2] For introductory material and references, see J.F. Gunion, H.E. 
Haber, G.L. Kane and S. Dawson, The Higgs Hunters Guide, 
Addison-Wesley Publishing. 



5 



[3] Phenomenological considerations for Higgs triplet representa- 
tions, especially in left-right symmetric models, are outlined in 
J.F. Gunion, J. Grifols, B. Kayser, A. Mendez, and F. Olness, 
Phys. Rev. D40 (1989) 1546; N.G. Deshpande, J.F. Gunion, B. 
Kayser, and F Olness, Phys. Rev. D44 (1991) 837. 

[4] Phenomenology for Higgs triplet models, focusing on the p — 1 
model with one \Y\ — and one \Y\ = 2 triplet, is considered 
in: J.F Gunion, R. Vega, and J. Wudka, Phys. Rev. D42 (1990) 
1673. 

[5] Recent work includes: M. Bando, T. Sato and T. Takahashi, 
Phys. Rev. D52 (1995) 3076; B. Brahmachari, Phys. Rev. D52 
(1995) 1; B. Brahmachari and R.N. Mohapatra, Phys. Lett. B357 
(1995) 566; E. Ma, Phys. Rev. D51 (1995) 236; D.-G. Lee and 
R.N. Mohapatra, Phys. Rev. D52 (1995) 4125. 

[6] K. Huitu and J. Maalampi, Phys. Lett. B344 (1995) 217. 

[7] K. Huitu, J. Maalampi, A. Pietila and M. Raidal, h ep- 
ph/9606311. 

[8] M. Lusignoli and S. Petrarca, Phys. Lett. B226 (1989) 397. 

[9] E. Accomando and S. Petrarca, Phys. Lett. B323 (1994) 212. 

[10] M. Kruse et at., FERMILAB-CONF-96/238-E (1996). 

[11] A.T. Goshaw et ai, FERMILAB-CONF-96/209-E (1996). 

[12] F Abe et al. , Phys. Rev. Lett. 76 (1996) 3070. 

[13] J. Ohnemus, Phys. Rev. D44 (1991) 1403; ibid. (1991) 3477; J. 
Ohnemus and J.F Owens, Phys. Rev. D43 (1991) 3626. 

[14] Particle Data Group, Phys. Rev. D54 (1996) 1. 

[15] C. Loomis et al., FERMILAB-CONF-96/232-E (1996). 

[16] J.F Gunion, Int. J. Mod. Phys. All (1996) 1551. 



6 



